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ABSTRACT 

It has been shown that the behaviour of primordial gas collapsing in a dark matter mini¬ 
halo can depend on the adopted choice of 3-body H 2 formation rate. The uncertainties 
in this rate span two orders of magnitude in the current literature, and so it remains 
a source of uncertainty in our knowledge of population III star formation. Here we in¬ 
vestigate how the amount of fragmentation in primordial gas depends on the adopted 
3-body rate. We present the results of calculations that follow the chemical and thermal 
evolution of primordial gas as it collapses in two dark matter minihalos. Our results 
on the effect of 3-body rate on the evolution until the first protostar forms agree well 
with previous studies. However, our modified version of GADGET-2 SPH also includes 
sink particles, which allows us to follow the initial evolution of the accretion disc that 
builds up on the centre of each halo, and capture the fragmentation in gas as well as its 
dependence on the adopted 3-body H 2 formation rate. We find that the fragmentation 
behaviour of the gas is only marginally effected by the choice of 3-body rate co-efficient, 
and that halo-to-halo differences are of equal importance in affecting the final mass 
distribution of stars. 
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1 INTRODUCTION 


According to the standard model of the primordial star for¬ 
mation process, the very first stars, the so-called Popula¬ 
tion III (or Pop. Ill) stars form within dark matter (DM) 
halos with virial temperature ~ 1000 K and ma sses ^10^ Mc^) 
that c ollapsed at redshift z ^ 25-30 or above (iHaiman et al l 
199fil: iTegmark et, al.1 Il997l: lAbel et,'^ I2OO2I: lOloved bonsl: 
BroinmfcLarsoiT 20041 : lBromm^t~am2009l ) . The hydrogen 


atoms in the DM halos combine with the small abundances 
of free electrons via H + e“ ^ H“+7, followed by H“ + H ^ 
H2 + e“, where the free electrons act as catalysts and are 
prese nt as residue from the epoch of recombination at z ~ 
1100 (jBarkana k, Loebll200ll : Oardi k, Ferrarall2005l k 


The typical fractional abundances of H2 of ~ I0~^ are 
sufficient to permit the primordial gas in the minihalos to coo l 
as it collapses (see, e.g. lSusa et al.flQ^ : lYoshida et al.ll2003l i 
via H2 rotational and vibrational line emis sion. This occurs 
until the gas reaches a temperature ~ 200K (jOrnukai Nishil 


llQQSl i at a density of ~ I0^cm“^, at which point the H2 en¬ 
ergy levels move into the local thermodynamical equilibrium 
(LTE) with the kinetic temperature of the gas, and the re¬ 
sulting cooling time becomes longer than the free-fall of the 
gas. The gas now heats up slightly as the collapse proceeds. 
This transition from cooling to heating with increasing den¬ 
sity sets a characteristic Jeans length, allowing the gas to 
fragment with Jeans mass of Mj ~ lOOOM© at typical tem- 

J >eratures of T ^ 200 K and number densities n ^ 10^ cm“^ 
Abel et alJUo^. l2002l: iBromm et aDll999l. l2002l) . 

Once the gas reaches a number density of ~ 10®cm“^, 
further H2 forma tion is possible via the 3-body reactions 
(jPalla et al.lll983l k 

H + H + H^H2 + H (I) 

H + H + H2 ^H 2 +H 2. (2) 

Unfortunately, the rate coefficients for the above set of re¬ 
actions are very uncertain in the temperature regime which 
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is applicable to Pop. Ill star formation, with the current es- 
tima tes in the literature s panning two orders of magnitude 
(see over Abell I2OO8I : iGlover SavinI l2009l for discus¬ 
sion). However, independent of the choice of rate coefficient, 
calculations of the gravitational collapse of primordial gas 
show that by a number density of ~ 10^cm~^, almost all 
the atomic hydro gen is converted to H 2 (|Yoshida et ahlliooil : 
iTurk et al.l[2009l b 


The sudden formation of H 2 has two interesting effects. 
First, the process is accompanied by heating, since 4.4 eV is 
released for every H 2 molecule that is formed. If the formation 
rate is fast enough, this heating can momentarily stall the col¬ 
lapse against collapse. On the other hand, the ability of the 
gas to cool is proportional to (among other things) the frac¬ 
tional abundance of H 2 , so the more H 2 is formed, the more 
the gas can cool. In practise, the interplay between these two 
phenomena is complex, and de pends crucially on the adopted 
rates for the 3-body reactions ([Turk et al.ll201lh . Under cer¬ 
tain conditions, a chemo-thermal insta bility can also occur 
that may promo te the gas to fragment (|Yoshida et ahllioOTl : 
iGreif et ffi1l2013l b 


In addition to any fragmentation induced by the chemo- 
thermal instability, recent simulations conclude that the disk 
around the first primordial protostar becomes gravitation¬ 
ally unstable and forms a multiple syst em with low-mass 


protostars instead of a single protostar ( 

Glark et al.1 

2008, 

201 Ibl: IStacv et al. 

l2010l:lGreif et al.ll201ll: 

Smith et al.1 

2011; 

Dopcke et al. 12013 

). This can occur on timescales as short 


as 10 to 100 years after the formation of the first object. 

I Glark et al.l (|2011bl j show that the fragmentation of this disc 
results from the fact that the disc is unable to accrete onto the 
protostar faster than it accretes from the in-falling envelope. 

The study by I Turk et ahl (|201l| j used both SPH and 
AMR simulations to conclude that the choice of 3-body-rate 
coefficients can introduce significant uncertainty into the ra¬ 
dial velocity, temperature and accretion rate, and therefore 
can affect the rate at which t he protostella r disc is fed. This 
has recently been verified bv iBovino et al.1 (|2014| j. who used 
bo th old r ates and newly calculated improved rate given 
by iForrevI (|2013l b to reach similar conclusion as I Turk et al.1 
(|201lh . Although they found compressed spiral arms, with 
strong density contrast between the arms and embedding 
medium, these simulations were unable to follow the evolu¬ 
tion beyond the formation of the protostar, and so could not 
address the long-term behaviour of the disc structure. 


However the degree to which the disc fragmentation de¬ 
pends on the chemical uncertainties in the 3-body H 2 for¬ 
mation rate has never been systematically tested. Previous 
studies did not include sink particles, and could not study 
the fragmentation that occurs once the first object is formed. 
In this work we focus on the effect that the choice of 3-body 
reaction rate has on the fragmentation of the primordial gas. 


2 NUMERICAL METHODOLOGY 


2.1 Simulation setup 


We investigate two minihaloshalos obtained from the cosmo¬ 
logical simulations of iGreif et al.1 (|201lL[201^j. which use d the 
hydrodynamic moving mesh code Arevo (ISoringelllioiOl b We 
use the snapshots from these simulations at a point when the 
central number density is just below 10® cm“® - comfort¬ 
ably before the onset of the 3-body reactions - as the start¬ 
ing point for this project. Interpreting the mesh-generating 
points of Arepo as Lagrangian fluids particles, we then use 
th e Arepo output as the initial conditions for our GADGET- 
2 (|SDringe]|l2005l j implementation. The modifications to the 
standard GADGET-2 include a time-dependent chemical net¬ 
work for primordial gas, and a sink particl es to capture the 
formation of collap sing protostellar cores ([Bate et al.1 1 19951 : 
IjapDsen et al.|[200a b A f uller descrip t ion of the code and its 
features can be found in iGlark et al.1 (|2011al li3j. The conver¬ 
sion from the moving mesh to smoothed-particle hydrody- 
namics (SPH) form alism is the same as that performed by 
ISmith et al.1 (|201ll j. and indeed the minihalos are identical, 
both in terms of the baryonic and dark matter content. 

The origin al Arepo simulati ons resolve the Jeans length 
with 128 cells (|Greif et al.1 12 OIli b at which point further re¬ 
finement is deactivated, resulting in roughly constant mass 
particles within central 1000 AU. However, to investigate 
the effect of 3-body H 2 reactions, here we use the intermediate 
Arepo snapshots (when the central number density n ^ 10® 
cm“®) as the starting point. Although the initial size of the 
minihalos are ~ 3 pc, all the fragmentation and accretion take 
place in the central region of the halos where the SPH particle 
masses are roughly 10~^ Mq. In GADGET- 2 simulation the 
mass resolution for 100 SPH particles (e.g. iBate BurkertI 
Il997l^ is pa M©. 

The chemical and ther modynamical trea tment of the gas 
in the Arepo simulations of iGreif et al.1 (|201l| j is also identical 
to those in our GADGET-2 implementation for the range of 
densities that are studied in this paper. More details of the 
initial conditions for this study are given below. 

Of the various rates available in the literature for the 
3-body H 2 formation reactions, we choose the extreme rate 
coefficients: the slowest of whi ch is lAbel et al. (|200^ (here¬ 
after ABN02) alongside that of lElower Harris (| 200 Ti l which 
provides the fastest of the rates (hereafter EH07). We have 
used the c ollisional-dis s ociati on rate in the LTE limit as de¬ 
scribed in I Turk et ^ (|201lh . Eor our SPH calculation we 
used a temperature-dependent value calculated via the prin- 
ciple of det ailed balanced. Eor the ABN02 runs, data from 
lOrell (|l987l j are used for low density (and low temperature, 
T < 300 K) cases, and a power law (described in Table 1) for 
high density and temperatures . Eor EH07 runs, we used the 
rates given in lElower Harris (|2007l b Eor the radiative cool¬ 
ing rate, we use Sobolev approximation to calcul ate it in the 
optical ly thi ck region, as describ ed in detail in iGlark et al.1 
(|2011al j and lYoshida et al.1 (|2006l j . Eull details of these rates 
are given in Table 1. 
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Ref 

3-body H 2 formation rates, A:3bh2 
(cm® s“^) 

Temp range 

ABN02 

ABN02 

FH07 

1.3 X 10-32(T/300)-°'38 

1.3 X 10-32 ©/300©i °o 

1.44 X 10-26©1-54 

(T < 300A) 
(T > SOOA) > 

H 2 collisional dissociation rates (cm^ s ^) 

ABN02 

1.0670825X10“^° x(T/ 11605)2-°^2 
exp(4.463/T/l 1605) x( 1 + 0 .2472 T/ 1 1605)4-512 


FH07 

1.38 X 10-4 ^ r-i °25exp(-52000/r) 



Table 1. Summary of 3-body H 2 formation and dissociation rates 
adopted here. 


2.2 Initial Conditions 

The simulations with halol start with a maximum cen¬ 
tral cloud number density (n) 10®cm“^ and minimum n of 
71cm“^. The initial temperature of the atomic hydrogen is 
470K (max) and 60K (min). Halol contains lOSOM© of gas 
and 690855 SPH particles. The maximum and minimum mass 
of SPH particles are 0.16 M© and 1.3 xl0“^ M©. This means 
that the numerical resolution of our simulations for halol for 
100 SPH particles is 1.3 xl0“^ M©. The simulations with 
halo2 starts with a maximum n of 10®cm“^ and minimum 
n of 85cm“^. The initial temperature of the atomic hydro¬ 
gen is 436K (max) and 54K (min). Halo2 contains 1093M© 
of gas and 628773 SPH particles. The maximum and mini¬ 
mum mass of SPH particles are 0.05M© and 1.4 xl0“^ M©. 
This means that the numerical resolution of our simulations 
for halo2 for 100 SPH particles is 1.4 xl0“^ M©. With al¬ 
together 6 million particles in the simulation, we resolve the 
halo gas and its velocity structure signifi cantly better tha n 
in the SPH models of our previous study (|Turk et al.ll201lh . 
However, we note that none of the existing numerical calcu¬ 
lations are able to fully resolve the turbulent cascade in the 
halo gas and do justice to the large Reynolds numbers (of 
10^ or abov e) that characterise the accretion flow (see, e.g. 
chapter 4 in iKlessen fc Glover!l20id ). 

Each of the minihalos are simulated with the two con¬ 
trasting 3-body reaction rates, resulting in 4 simulations in 
total. The different minihalos allow us to distinguish between 
those features of the high-density collapse/evolution that are 
caused by the chemo-thermal differences arising from the re¬ 
action rates, and those features which result from cosmic vari¬ 
ance (i.e. differences caused by the formation of the miniha¬ 
los). 

As collapse progresses in the central region of the cloud, 
sink particles are created once the number density of the gas 
reaches 5 xl0^^cm“^, at which point the gas has a temper¬ 
ature of around lOOOK. The corresponding Jeans mass at 
this density and temperature is 0.06M©, so our calculations 
for both of the halos are well resolved. We replace a can¬ 
didate particle by a s i nk particle (see, e.g . iBate et ^Il995l : 
iKrumholz et al.l [20041 : Ijaposen et al.l[2005 i that can accrete 



Figure 1. Radially binned, mass-weighted averages of physical 
quantities for different 3-body H 2 formation rate (red: ABN02, 
green: FH07) of two different cosmological halos (dotted: halol, 
solid: halo2), are compared just before the first protostar forms. 
The enclosed mass (A), H 2 fraction (B), temperature (C), ratio 
of cooling time to free-fall time (D) and ratio of free-fall time to 
sound-crossing time (E) are plotted as a function of number density. 

gas particles within its accretion radius race that we fix at 6 
AU, which is the Jeans radius at the density threshold for sink 
creation. We also prevent sink particles from forming within 
2 race of one another in order to avoid spurious formation of 
new sink particles out of gas that, in reality, would by ac¬ 
creted by a neighbouring sink particle. Lastly, gravitational 
interactions between sinks, and between the sinks and the gas, 
are softened using a fixed softening parameter of 1.2 AU for 
the sinks. We run our simulations until a mass of 21 M© has 
been accreted by all the sink particles.The final time scales 
range between 2500-4000 yr due to the physical differences 
between the minihalos and the 3-body H 2 formation rates. 


3 PHYSICAL CONDITIONS JUST BEFORE 

SINK PARTICLE CREATION 

In this section we examine the physical conditions in the gas 
once the central region has collapsed to a density of ~ 5 x 
10^^ cm“^. The radial profiles of the gas for the two different 
halos and two different 3-body rates are shown in Figure [1] 
The panels show mass-weighted averages of the properties of 
individual SPH particles within radial logarithmic bins. 

Panel A shows the mass enclosed as a function of the 
number density, and we see that it displays a rough power- 
law over most of the density range covered here. This simply 
follows from the radial power-law in the de nsity which is typ- 
ical for primordial gas: n{r) oc (e.g. lAbel et al.ll2002l: 
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lYoshida et, al.l 1200(11 . l2008h . We see that there are slight dif- 
ferences between the profiles for the two rates, but these are 
also comparable to the differences between the two halos. 

In panel B, we see the rapid rise in the H 2 fraction with 
increasing density that arises from the 3-body reactions. Here 
we see a clear difference between the two rates, with the runs 
with FH07 producing most of the H 2 by a density of 10^ 
cm“^, while the ABN02 runs are only starting to produce 
H 2 by the time this density is reached. As the rates are a 
strong function of density, we see very little difference in the 
H 2 fractions between the two halos. 

The mean temperature as a function of the number den¬ 
sity is shown in panel C. We find that the gas is significantly 
cooler, up to a factor of 2, around density ~ 10^^ cm“^, in 
the FH07 halos than in those using ABN02. The differences 
in temperature-density evolution start around a number den¬ 
sity of 10^ cm“^, once the H 2 abundance in the FH07 case 
has increased by an order of magnitude. Again, we find that 
there is little difference between evolution of halos 1 and 2 
in this panel, implying that the chemical state of the gas is 
more important for the temperature evolution than any dif¬ 
ferences in the halo properties. C urves in panels A , B and C 
are consistent with the results of iTurk et al.l (|201lh . 


4 VELOCITY STRUCTURE AND ACCRETION 


In this section we study the difference in mass accretion as¬ 
sociated with the cloud collapse, that can arise due to the 
choice of 3-body reaction rates. In what follows, we will first 
study the accretion rate as predicted by the properties of the 
gas in the clouds just before the first sink particle forms. We 
will then examine the accretion rate as measured by the sink 
particles themselves. 

The accretion rate can be estimated from the radial gas 
profiles using, 

M(r) = 4 7rr^p(r)t;rad(r). (3) 


This has been shown to be a good estimate of the accretion 
rate (|Clark et al.l 12011^ 1. In addition, one can define an ac¬ 
cretion time from. 


Tfenc ) 
47rpi;rad ^2 ■ 


(4) 


Since accretion is related to infall velocity, we also discuss 
the corresponding radial and rotational velocities. Figure [2] 
shows these physical quantities as the function of enclosed 
mass for all the simulations, using the radial profiles in the 
gas just before the formation of the first sink particle. The 
left panel of Figure [2] in the clockwise order, shows the ac¬ 
cretion rate, radial velocity, radial velocity over sound speed 
and accretion time respectively. Similarly the right panel of 
Figure [2] shows the rotational velocity, radial velocity over ro¬ 
tational velocity, rotational velocity over Keplerian velocity 
and rotational velocity over sound speed respectively in the 
clockwise order. 

We note that the mass accretion rate (M) for all the 
simulations has a maximum at an enclosed mass ~ 1-2 Mq. 


Given that H 2 line-cooling becomes optically thick at the cor¬ 
responding densities for this enclosed mass, it makes sense 
that M for all the simulations converge at this mass scale, 
as the gas loses its ability to cool efficiently. In general the 
differences in accretion rates caused by the choice of 3-body 
formation rate are small, around 25 percent. The differences 
between the rates in different halos are much larger, at around 
a factor of 2 for masses below a few tens of solar mass. How¬ 
ever we see from the top left panel of Figure [2] that once 
the mass approaches 100 M©, the difference in the accretion 
rates can be large, approaching an order of magnitude. We 
conclude that halo-to-halo variations will play a larger role in 
the accretion rates than chemical uncertainties. 

The choice of 3-body rate coefficient introduces differ¬ 
ence in radial velocity of ~ 0.4-0.5 km . We also note 
that the radial velocities for both rates are substantially lower 
than the rotational velocities. The gas particles spiral inwards 
with trajectories following a combination of gravitational at¬ 
traction and angular momentum conservation. All clouds are 
sub-Keplerian. However, there is a slight tendency for the gas 
to be more rotationally supported in the FH07 runs - at least 
for the inner 10 M© of gas. 

We can compar e these resu l ts wit h both the AMR and 
SPH simulations of iTurk et al.l (|201lh . Our high-resolution 
study also concludes that the difference in the 3-body H 2 
formation rates significantly affect the dynamical behaviour 
of the collapsing gas. The uncert ainty in the velocity struc¬ 
ture has also been investigated bv iBovino et al.l (|2014h using 
hydrodynamics code ENZO. They also found that there are 
relatively large differences in the radial velocities of about 2-3 
km s~^ depending on the halo distribution as well as 3-body 
H 2 formation rates. 

Although there is scatter in the velocities for ABN02 and 
FH07 rates, the scatter also exists in the velocities for both 
halol and halo2. We again conclude that the halo-to-halo dif¬ 
ferences in the velocity structure is, in general, more pro¬ 
nounced than that caused by the differences in the chemical 
reaction rates. 


5 FRAGMENTATION 

We have shown in Section [22] that the temperature of the gas 
at high densities, such as those in disc, depends on the choice 
of adopted 3-body H 2 formation rate. In this section, we study 
how the gas fragments as it evolves beyond the formation of 
the first collapsing core. Our discussion will first focus on 
the properties of the gas, to see whether there is any hint 
of the future evolution already present before the first sink 
particle forms. We will then compare the properties of the 
sink particles in the simulations and discuss the differences. 

One useful measure of the ability of the gas to fragment 
is the ratio of the thermal to gravitational energy. This can 
be represented by the number of Bonnor -Ebert masses (Mbe) 
conta ined in the central dense volume (|Eber i] ll95,4: lB onnorl 
Il956[ ). We therefore investigate the changes in the number 
(i.e, Menc/MBE) for each of our simulations and how it varies 
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Figure 2. Radial logarithmic binned, mass-weighted averages of mass accretion rate (A), radial velocities (B), ratio of radial velocity 
to sound speed (C), accretion time (D), rotational velocity (E), ratio of radial to rotational velocity (F), ratio of rotational to Keplerian 
velocity (G), and ratio rotational velocity to sound speed (H) are plotted as a function of enclosed mass. The halo to halo variation as well 
as the uncertainty in the 3-body rate introduce an uncertainty in the mass accretion rate, radial and rotational velocity. 



Figure 3. The number of Bonnor-Ebert masses (Menc/^BE) in 
the central dense volume is plotted as a function of enclosed mass, 
just before formation of first protostar. The degree of gravitational 
instability for FH07 within the central dense regime is larger than 
ABN02, indicating larger opportunity for the gas to fragment in 
the FH07 case, than in the ABN02 simulations. 


wi th the enclo s ed ga s mass, similar to the analysis performed 
bv lAbel et al.l (|2002l j . We computed the Bonnor-Ebert mass 
as the mass-weighted average within the logarithmic radial 
bin, 

Mbs = (5) 

where Cs is the sound speed, Pext is the external pressure 
that we assume to be equal to the local gas pressure, /x is the 


mean molecular weight and 7 = 5/3 is the adiabatic index 
respectively. 

The results are shown in Fig. O Note that the ra¬ 
tio Menc/MsE is equivalent to the ratio of fragmentation 
timescale to that of accretion, tuag/^acc- Once again, these 
snapshots were taken immediately before the formation of 
the first protostar. The dashed lines represent the case when 
fragmentation timescale equals the accretion timescale. 

While we see that there are differences between the mass 
profiles of the FH07 and ABN02 runs, the effects are, once 
again, not substantial. However, in the case of Menc/MeE 
we see that the FH07 runs have a differently shaped profile 
from the ABN02 runs (top panel of Fig. [3]). fn the FH07 case, 
the profiles are peaked at characteristic mass of 1-2 M©. In 
contrast, the ABN02 runs do not display such a peak, and 
it is less clear whether these runs would favour a particular 
mass for fragmentation. 

If the ratio Menc/MeE < 1 the gas enclosed in the shells 
are accreted faster than it can fragment. As a result fewer 
new protostars are formed and the avail able mass contribute s 
to the mass growth of the existing ones (|Dopcke et al.|[201^ i. 
In our case we find Menc/MeE > 1, i.e., the gas in the shells 
can fragment faster than it is accreted by the central dense 
clump, favouring low-mass protostars. 

Therefore by considering the radial density profiles before 
the onset of sink formation we conclude that gas in FH07 case 
is perhaps more likely to fragment than in the ABN02 case. 
Next we study how the gas fragments at later epochs after 
formation of the first sink particle. We follow the simulations 
until 21 M 0 have been converted into, or accreted onto, sink 
particles. 

We have calculated the total mass accretion rate by 
all the sinks particles. Although fragmentation may depend 
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Figure 4. Surface density and temperature in a region of 2000 AU centred around first protostar for different 3-body H 2 formation rates, 
are shown when the total mass of the sinks reaches ~ 21 Mq. 


on the way sink particles are modelled, our implementation 
(Section 2.2) takes considerable care to avoid what is of¬ 
ten called artificial fragmentation and the spu rious forma¬ 
tion of sink parti cles (further deta il s are found in Clark et al.l 


201 Ibl. see also iBate et al.l Il995l . Ijaposen et al 
Federrath et al.l[2010ll. 


Clark e 

l2005l. 


and 


The column density and column-weighted temperature 
distribution in the inner 2000 AU at the end of the simula¬ 
tions are shown in Figure [H In all cases we see that the sim¬ 
ulations exhibit disc structure on several scales within this 
central region. It is not surprising, given the high levels of 
rotational support seen in right-hand panel in Figure O We 
also see that the temperatures in the disc are slightly less in 
the case with the FH07 rate, due to the increased H 2 fraction. 
However from these images it is clear that ah simulations frag¬ 
ment to form a small-A" system within the time that 21 Mq 
of material is accreted onto sink particles. 

Figure [5lA shows the time evolution of total mass accre¬ 
tion rate onto the sinks (same as Figure [21). The plot shows 
that mass accretion rate for ABN02 is larger than the FH07 
rate. For both rates, dM*/dt decreases with time. After 500 
years, further sink particles form, and the total accretion rate 
increases again, however, now with large temporal variations. 
In Figure [5p, we plot the mass accretion rate as a function 
of total sink mass. Figure [5p shows the time evolution of the 
total mass in sinks and Figure [5p shows the maximum sink 
mass in the period over which the sink formation occurs. 

We see that the simulations employing the ABN02 rate 
produce a maximum sink mass that is consistently above the 
FH07 rate, for a given halo. However, the difference is less 
than a factor of 2 and is similar to the inter-halo difference. 


When we examine Figure [5fl, which shows the total number 
of sink particle versus the total mass in ah sinks, we see a 
similar behaviour: the differences between the runs with dif¬ 
ferent rates can be a factor of two. Although in this case the 
halo-to-halo differences are less pronounced. 

We caution the reader that our simulations do not in¬ 
clude feedback effects from nascent protostars. We therefore 
restrain ourselves from following the dynamical evolution of 
the embedded stellar system over very long timescales. We 
end our calculations before we expect that radiative and me¬ 
chanical feedback changes the thermal as well as chemical 
state - and as a consequence - also modihes the fragmenta- 
tion behavior of the remaining gas (see , e.g., the discussion in 

IWise AbefcOOSl.IWhalen et aLllioi^ . iHosokawa et al.|[201lL 

or lStacv et ^|20ld). 


6 SUMMARY 

We have investigated the effects of different proposed 3-body 
H 2 formation rates during the collapse of primordial star¬ 
forming clouds and analysed their inhuence on the result¬ 
ing fragmentation behaviour of th e gas. We compared the 
rates proposed bv lAbel et al.l (|2002l l and bv [Flower Har^ 
(|2007f ). which span two orders of magnitude. We follow the 
chemical and thermal evolution of the primordial gas in two 
different minihalos in order to assess differences introduced 
by varying halo parameters. We hnd that the uncertainty in 
the 3-body H 2 formation rates leads to the differences in the 
chemical state of the gas which is more important for the 
thermal evolution of the clouds than any differences in the 
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halo properties. The halo to halo variation affects the dy¬ 
namical evolution of the gas and mass-accretion more than 
the chemical uncertainties. 

We a l so fin d that the simulations employing the 
lAbel et H] (|2002l l rate to produce on average fewer and more 
massive fragmen t s com pared to those calculations using the 
Flower Harrisl (l2Q07l l rate (consistent with the results of 
Bovino et al .1 il2ni4ir which however did not use sink parti¬ 
cles). This difference of roughly a factor of two, however, is 
comparable to the hal o-to-halo variation . This is similar to 
the results reported bv lTurk et al.l (|201lh with added advan¬ 
tage of using sink particles. 

In summary, we find that the variations introduced to 
the mass distribution of primordial stars introduced by the 
uncertainties in the three-body H 2 formation rate is of similar 
order to the fluctuations introduced by difference in the halo 
properties. 
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Research Gouncil under the European Gommunity’s Seventh 
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Figure 5. Time evolution of the protostellar system: Accretion rate as the function of time (A) and total mass in sink particles (B). The 
total mass of all the sinks particles (C) and the most massive sink particle (D) are plotted as function of time. The simulations employing 
FH07 produce more fragmentation than ABN02 (E). 
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